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Abstract 

A general exact analytic expression for the probability of matter-enhanced 
two-neutrino transitions in a medium (MSW, RSFP, generated by neutrino 
FCNC interactions, etc.) is derived. The probability is expressed in terms 

£> | of three real functions of the parameters of the transitions: the "jump" 

-^ \ probability and two phases (angles). The results obtained can be utilized, 

in particular, in the studies of the matter-enhanced transitions/conversions 
of solar and supernova neutrinos. An interesting similarity between the 
Schroedinger equation for the radial part of the non-relativistic wave function 
of the hydrogen atom and the equation governing the MSW transitions of 

Q_i| solar neutrinos in the exponentially varying matter density in the Sun is also 

Q_i' briefly discussed. 

> 
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1. Introduction 

In the present article we derive exact, simple and general analytic expression for the 
probability of matter-enhanced two-neutrino transitions in a medium. The matter-enhanced 
neutrino transitions provide, as is well-known, at least three different varieties of neutrino 
physics solutions of the solar neutrino problem Jl|-|5[. Such transitions, in particular, can 
play important role in the supernovae dynamics 0J^] and can be at the origin of the ob- 
served large space velocities of pulsars ||. Even in the simplest case of transitions involving 
only two weak-eigenstate neutrinos, however, the system of evolution equations describing 
the transitions does not admit, in general, exact solutions and one has to rely on numeri- 
cal methods to calculate the corresponding transition probabilities. There are few notable 
exceptions: the evolution equations can be solved exactly, for instance, in the case of MSW 
transitions |l],|2| in matter with density which changes linearly [p|-|TTH (see also fl2|| ) or ex- 
ponentially [0,0] along the neutrino path. On the basis of the exact solutions expressed 



in the two cases respectively in terms of parabolic cylinder and confluent hypergeometric 
functions, simple expressions for the neutrino transition probabilities, containing only ele- 
mentary functions, have been obtained [] [p|- [TT|JT3|] . The case of exponentially varying density 
is especially relevant for the analytic description of the solar neutrino transitions in the Sun 
since according to the contemporary solar models JI7| the density decreases approximately 



exponentially from the center to the surface of the Sun. Indeed, it was found that the expres- 
sion for the average probability obtained in the exponential density approximation provides 
a very precise (and actually, the most precise) description of the MSW transitions of the 
solar neutrinos in the Sun ||TH| . 



The expression for the average probability derived in the linear density approximation 
contains as an integral part the result by Landau and Zener for the corresponding "jump" 
probability fL2l]16| . It is widely used, for example, in the studies of the MSW transitions 



and resonance spin-flavour precession (RSFP) or conversion || of the supernova neutri- 
nos (see, e.g., [§]-§])• However, its accuracy in describing the supernova neutrino transi- 
tions/conversions has never been tested. More specifically, the oscillating terms present 
in the relevant transition/conversion probabilities have always been neglected without any 
control on the validity of this approximation because no explicit expressions for these terms 
could be derived Q 

The analytical result for the two-neutrino transition/conversion probability we obtain in 
the present study, in addition of being exact and rather simple, is also general and universal in 
form. It is valid for two-neutrino MSW transitions, RSFP, for matter- enhanced transitions 



x The expression for the average MSW neutrino transition probability in the case of linearly varying 
density was derived first in ref. |15| using qualitative arguments and the Landau - Zener formula 
pUTl for the "jump" probability. 



2 One exception is the first article in ref. pi] where the results of the present study have been used 

to ensure that the corresponding oscillating terms are strongly suppressed. 
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induced by neutrino flavour- changing neutral current (FCNC) interactions @J^], etc. In 
all these cases the neutrino transition/conversion probabilities are shown to be given by 
expressions which have one and the same universal structure. Such an expression was derived 



first in ref. ]T9] for the solar neutrino MSW transition probability in the exponential density 



approximation. 

2. The Probability of Two-Neutrino Transitions in a Medium: General Results 

Consider neutrino transitions in a medium, which involve two weak-eigenstate neutrinos 

u a and up, u a — > up, u a ^ up — v e , u fl , u T , u s , v s being a sterile neutrino f\ We shall 
assume that the transitions are described by a system of evolution equations which has (or 
can be reduced to) the form: 

d/A*(Mo)\ f-e(t) e'(t)\ (A a (t,t )\ 
l dt \ Ap(t,t ) ) ~ { e'(t) e(t) ) { Ap(t,t ) ) [i) 

where A a (t, t ) (Ap(t, to)) is the amplitude of the probability to find neutrino v a [up) at time 
t of the evolution of the neutrino system if at time to the neutrino u a or vp has been produced, 
t > to, e(t) and e'(t) are real functions of some of the physical quantities characterizing the 
medium (density, magnetic field strength, etc.) and the neutrino system in vacuum (neutrino 
energy, mass squared difference, etc.), and we have omitted the indices a and (3 from e(t) and 
e'(t). In writing the evolution equations in the form (1) we have supposed that the effects 
of neutrino absorption, creation and of possible neutrino instability (or loss of coherence) 
in the evolution of the neutrino system are negligible. This corresponds to a large number 
of physically interesting cases, in particular, those discussed in [BD-H]. Under the above 
assumption the evolution matrix of the system is hermitian and the equations describing the 
two-neutrino transitions in a medium can always be brought (by phase transformations of 
the two probability amplitudes) to the form (1) with real e(t) and e'(i). 

There are two specific types of initial conditions for the system (1) relevant to our dis- 
cussion: 

A° a = l, A, = 0, (A) 

and 

< = 0, 4 = 1, (B) 

where A° a = A a (t ,to) and Aa = Ap(to,t ). If the initial conditions (A) hold, A a (t, t ) and 
Ap{t, t ) are probability amplitudes of the u a survival and of the v a — > up transition while the 
neutrinos propagate from the point of their production to the point of neutrino trajectory 
reached at time t: A a (t,t ) = A{y a — * v a ) and Ap(t,t ) = A{y a — > vp) (we have not 
indicated explicitly that A(u a — > v a ) and A{y a — > vp) depend on t and to)- Similarly, if the 
initial conditions (B) are valid we can write A a (t,to) = A{yp — > u a ), Ap(t,t ) = A{yp — > up), 



3 In what follows the term "transition" will be used to denote the process v a — > up of the change 
of the type of the weak-eigenstate neutrino for all possible different cases, including the case of 
RSFP. 



in accordance with the interpretation of the probability amplitudes A a (t, to) and Ap(t, to). 
In the case of most general initial conditions 

<^o, aj^o, i<r+i^i 2 = i, (2) 

the solutions of the system (1) are expressed as linear combinations of the solutions corre- 
sponding to the initial conditions (A) and (B): 

A*(t, t ) = A° a A{u a -► u a ) + A° p A{vp -► u a ), (3a) 

A (t,to) = A° a A(u a - v p ) + A° A(up -> u ). (36) 

Evidently, for any of the initial conditions (A), (B) or (2) one has in the case under study 

\A a (t,t )\ 2 + \A /3 (t,t )\ 2 = l J (4) 

which implies 

A(v a -► v a )A*(vp -* v a ) = -A*{vp -> vp)A{v p -> i/ a ) . (5) 

We give next several concrete examples of matter-enhanced neutrino transitions in a 
medium to which our general results will apply. In the case of MSW v e — > v^{ T ) transitions 
of solar or supernova neutrinos we have [|I|,pO| 

where Am 2 = m| — m 2 , m lt 2 being the masses of two neutrinos h> lt2 with definite mass in 
vacuum, E is the neutrino energy, 9 is the neutrino mixing angle in vacuum, and N e (t) is the 
electron number density at the point of neutrino trajectory in the Sun or supernova, reached 
at time t. For the transitions z/ M ( r ) — > v e which can take place in supernovae, e(t) and e'(t) 
are also given by (6). If the MSW transitions are into sterile neutrino, v e — > u s , N e (t) in eq. 
(6) has to be replaced by |H] (N e (t) — ^N n (t)), where N n (t) is the neutron number density. 
In supernovae MSW v^m —* v s transitions can also take place. In the latter case N e (t) in eq. 
(6) has to be substituted by (— |iV n (£)). Finally, for the corresponding transitions involving 
antineutrinos the term with the Fermi constant in e(t) changes sign. 

Another relevant example is the simplest version of neutrino resonance spin-flavour pre- 
in the Sun or in supernovae. It is described by (1) with || 

V2G F (N e (t) - N n (t))}, e'(t)=fi u B ± (t). (7) 

m 2 {p^ T )) — m 2 (z/ e ), fi u is a v e — v^) transition magnetic 
moment of a dipole type and B±(t) is the value of the component of the solar or supernova 
magnetic field perpendicular to the neutrino momentum. The RSFP of the type u e — > i^( T )? 
which was shown in ref. to be one of the possible mechanisms of supernova shock revival, 
is also described by (1) with e(t) and e'(t) given by eq. (7) in which the term containing the 
Fermi constant has an opposite sign. For neutrino RSFP in a twisting magnetic field [f2~2~| , f2~3| 
as well as MSW transitions or RSFP of supernova neutrinos, which can be responsible for 



cession, e.g., v e -> u^ T ), 


in the 


e(t) = - 

v ; 2 


r Am 2 
^ 2E 


In eq. (7) Am 2 = m 2 ■ 
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the observed large space velocities of pulsars ||, the expressions for e'(t) coincide with those 
given in eqs. (6) (MSW transitions) and (7) (RSFP), while the expressions (6) and (7) for 
e(t) are modified (but not essentially from the point of view of the problem we are interested 
in) by the presence of additional terms 0. 
Finally, the system (1) with 

e(t) = \V2G F N e {t) [A(l + 2 ML) - 1], e\t) = s/2G F N e (t) A'(l + 2 ^||), (8) 



where A and A' are real constants, A > 0, corresponds to matter-enhanced transitions of 

solar {y e — > u T ) or supernova ( v e — > v T ^, v r^) - * v e ) neutrinos induced by v e flavour- 
changing and new v e and v r U) flavour-conserving but flavour non-symmetric, neutral current 

interactions (on the d— quark) ||. Note that in the latter case matter-enhanced v e — > v T ^ 

( v t(u)~* v e) transitions are possible even if neutrinos have zero mass and neutrino mixing 
is absent in vacuum. 

Our aim is to derive an exact and general expression, e.g., for the v a survival probability 
assuming that the functions e(t) and e'(t) are given f\ This is typically the case in the 
studies of the effects of the matter-enhanced transitions of the solar and supernova neutrinos, 
indicated above. In these studies one uses the predictions for the electron and neutron 
number density distributions and for the neutrino energy spectra provided by the standard 
solar models and by the models of the supernovae. The energy spectrum of solar neutrinos 
is practically solar model independent and is known with a relatively high precision [ |T7|j . 
The supernova neutrino energy spectra are somewhat model dependent |J|7|]. No direct 
information exists about the magnetic fields in the solar and supernovae interiors, and in 
the studies of the neutrino conversion in which they are presumed to play important role, 
one uses plausible field configurations (see, e.g., the third article quoted in ref. || and ref. 
J7j). The magnetic field strength and Am 2 in the case of RSFP, the neutrino mixing angle 
in vacuum 9 and Am 2 in the case of MSW transitions, and the constants A and A' in the 
example specified by eq. (8), are often treated as free parameters to be determined by the 
physics of the problem being investigated. We will not specify the form of e(t) and e'(t) in 
our further analysis. Our results will be valid for any e(t) and e'(t) and, in particular, for 
e(t) and e'(t) given in eqs. (6) - (8) and the cases discussed above. 

Let us remind the reader that neutrinos go through a resonance point reached at time 
t res , to < t res < t, on their way to the final point of their trajectory if the function e(t) 
changes sign going through zero at t — t res , e{t res ) = 0, while e'(t res ) ^ (see the first 
article quoted in ref. |2(J as well as ref. |J| and, e.g., [|,H,^|). These conditions are typically 



fulfilled in most of the cases of physical interest, in which also one has e'(t) ^ on the whole 
neutrino trajectory, except possibly at the trajectory's final point. For v e -^ u^ T ) MSW 



4 In all the indicated cases the neutrinos are assumed to be relativistic. 

5 The results we shall obtain can be relevant also for the inverse problem, namely, the problem of 
reconstructing e(t) and/or e'(t) from data about v a (y a — > vp) survival (transition) probability. 



transitions, for instance, the resonance condition can be satisfied, as it follows from eq. (6), 
if, e.g., Am 2 > and COS26 1 > 0, and these two inequalities will be assumed to hold in our 
further discussion of the MSW effect for solar neutrinos. 

It is convenient to introduce the neutrino mixing angle in matter, 

e'(t) 
sm29 m (t)= , = (9) 

and the neutrino matter-eigenstates at time t, ^12 (t): 

|zy Q > = \u?(t) >cos9 m (t) + \i/?(t) > sin0 m (t), (10a) 

[vp > = -\v?(t) > sm9 m (t) + \v™(t) > cos9 m (t). (106) 

At the resonance point we have: sin 2 29 m (t res ) = 1. The states ^{^(t) >, which are also 
called "adiabatic" , are the instantaneous eigenstates of the evolution matrix (Hamiltonian) 
in (1) at time t, corresponding to the two eigenvalues, E™ 2 (t), whose difference is given by 



E™-E™ = 2 v /e 2 (£) + e' 2 (£). (11) 

Under certain specific conditions (e.g., u e produced in the Sun at densities N e (t ) much 
larger than the resonance density N™ 3 = Am 2 cos 29/(2 E\t2Gf)) the states (^[^(to) >, for 
example, can practically coincide with the states \vp, a >■ 

We shall denote by A(v|™(i ) — ► vT(t)) = A^(t,to) the probability amplitudes of the 
v T{to) ~~ * u T(t) transitions, i,j = 1,2, which take place when the neutrinos propagate in 
the medium. The "jump" (or "level crossing") probability, P', 

P' = \A(y?(to) - K(t))\ 2 = |A-| 2 , (12) 

plays a very important role in the neutrino matter-enhanced transitions in a medium: its 
value determines the type of the v a {8) ~^ u /3(a) transition (for P' = it is adiabatic and 
nonadiabatic otherwise) and typically controls the value of the transition probability in a 
large region of the corresponding parameter space Q. Moreover, in a wide class of cases one 
can use the Landau-Zener result for P' [|l|,|16|] or, e.g., its analog derived in the exponential 



density approximation [|J]. Then P' is determined by the values of e(t), e'(t) and of their 
derivatives at the resonance point (see further). Therefore we would like to find an exact 
general expression for the v a {p) survival and v a (p) ~^ u i3(a) transition probabilities in terms of 
the "jump" probability. To this end let us express the four probability amplitudes A(v a — > 
v a (j3)) an d A(v$ — > v a (p)) i n terms of the amplitudes A™, i,j = 1, 2. This can be done using 
the relations (10a) and (10b): 

A\y a > v a ) = A^c m c m + /i 12 c m s m + /i 21 s m c m + /i 22 s m s m , [loa) 



6 As we shall see later, one has in the case of interest: \A(u2 l (to) — > ^(t))] 2 = {Aty™ (to) 
v?(t))\ 2 =P>. 



A[y a ► up) = —A-^c m s m + A 12 c m c m — A 21 s m s m + A 22 s m c m , {L6b) 

Ayvp — > v a ) — ~A 1 iS m c m — A 12 s m s m + A 21 c m c m + A 22 c m s m , (13c) 

A\Vfi — * ^/3j = ^-ll S m S m ~~ J ^12 S m C m — A 21 C m S m + ^l22 C m, C m, (13a) 

where c$ = cos6> m (£(o)) and s^ = sin0 m (£( O )). 

For any fixed to and t > to one has from (4) and (13a) - (13d): 

I A m \ 2 -l- I A m \ 2 — 1 I 4 m l 2 -l- I A m \ 2 — 1 A m ( A m Y — —( A m Y A m (~\A\ 

\ A XX\ + 1^121 — X > l /1 2ll + l /1 22l — l i /1 lll /1 2lJ — 1^22J ^12 ■ l i4 J 

We shall prove next that the solutions of the system of equations (1) corresponding to 
the initial conditions (A) and (B) satisfy the following relations: 

A*{y a -> v a ) = A(v0 -> z^), A(^ -» i/ a ) = -A*(z/ a -> up) . (15) 

Indeed [], using the first equation in (1) to express Ap{t, t ) as 

A {t,t o ) = ~ (e + i-)A a (t,t ), (16) 

and substituting eq. (16) in the second equation in (1) we get a second order differential 
equation for A a (t,t ): 

where e^ = -^^'\ It is easy to show that A*Jt,to) satisfies exactly the same second order 
differential equation: this can be checked, e.g., by taking the complex conjugate of the two 
equations in (1) and by using the second equation written as 

-A* a (t,t ) = ^ (e + ij t )A}(t,t ), (18) 

to eliminate (— A* a {t, to)) from the first equation. Moreover, the initial conditions (A) lead 
to initial conditions for the solution A a (t, to) of the equation (17), which coincide with the 
initial conditions for the solution A*^(t,t ), following from conditions (B). This fact and the 
relations (16) and (18) lead to eq. (15). 

It is easy to convince oneself utilizing eqs. (13a) - (13d) and (15) that the amplitudes 
A™ satisfy analogous relations: 

A™ = (A? i y, A™=-(A? 2 y. (19) 



7 For the special case of MSW v e — v kV T ) transitions, eq. (6), these relations were shown to be 
valid in refs. HH). 

8 Obviously, the proof of eq. (15) presented here relies on the assumption that e'(t) ^ and that 
the derivatives of e(t) and e'(t) exist on the neutrino trajectory. 
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These relations can be derived also from the systems of evolution equations for the amplitudes 
Aii, ^12 an d f° r the amplitudes A™ x , A™ 2 which can be obtained from (1) by using (13a) - 
(13d). It follows from (19), in particular, that: 

$11 (Mo) =" $22 (Mo), $2l(Mo)=7T-$ 12 (Mo), (20) 

where $y(Mo) is the phase of the amplitude A^{t,t ), 

arg(A%) = $ ij ,i,j = 1,2. (21) 

Thus, eqs. (14) and (19) imply that the four complex amplitudes A™, i,j = 1,2, depend 
actually only on three real functions of the parameters of the problem, which can be chosen 
to be the "jump" probability P' = \A^\ 2 = |v4^|| 2 and the phase functions $12 (Mo) an d 

$22 (Mo)- 

It is not difficult to find now the expression of interest for the v a survival probability, 
P(y a — ► u a ), assuming the initial conditions (A) are valid: P{y a — > u a ) = \A(u a — > Vq)\ 2 . 
Note that due to eqs. (4) and (15) the other three relevant probabilities, P{yp —* vp) = 
\A(i/0 — > z^g)| 2 and P{y a (p) — > vp(a)) — \A(v a (p) ~~ > ^(a))! 2 , can all be expressed in terms of 
P{y a — » u a ). Using eqs. (12), (14), and (19) - (21) we get from eqs. (13a): 

4 
P{y a -> v a ; t, t ) = P{u a -► v a ) + Y, P r° SC (t> t ), (22) 



r=l 



where 



p(y a _ Va ) = I + (I - p>) cos29 m (t ) cos29 m (t), (23) 



is the average probability and 



Pr(Mo) = V^ 1 - P ') c °s2^ m (to) sin2# m (t) cos($ 12 + $ 22 ), (24a) 



P° sc (t,t ) = -^P'(1-P')sm26 m (t )cos26 m (t) cos($ 12 - $ 22 ), (246) 

P 3 ° sc (Mo) = -~P'sin20 m (t o ) sin2^ m (t) (cos2$ 12 + cos2$ 22 ), (24c) 

^T^ *o) = 2 sin 2e ^o) sin 20m(*) cos 2$ 22 , (24d) 

are oscillating terms. The functions $i 2 (t, to) an d $22(^5 *o) are responsible for the oscillatory 
dependence of the four probabilities P(v a (3) — > v a (p)) an d P(. u a((3) ~^ u /3(a)) on the neutrino 
energy E and/or on the parameters characterizing the transitions being studied. 

The result expressed by eqs. (22) - (24d) was derived for the probability of MSW transi- 
tions (see eq. (6)) of solar neutrinos in the Sun in ref. [[HJ, utilizing the exact solutions of the 
system (1) found in the case of electron number density N e (t) changing exponentially along 



the neutrino trajectory f\. Explicit analytic expressions for the "jump" probability P' as well 

as for the phase functions $i2(Mo) an d $22 (Mo) m the indicated case were also derived 

T3| , |T9|| (see also f26||). This permitted to perform a detailed study of the magnitude and 



the behavior under various averagings of the oscillating terms present in the solar neutrino 
MSW transition probability (for details see refs. |19],|27|) Q. 

For adiabatic transitions one has P' = and consequently P°| c 3 (t,to) — 0, while 
P£ sc (t,to) can be non- negligible. Thus, -Pff^Mo) can be identified as nonadiabatic os- 
cillating terms, while for P' = 0, P£ sc (t,to) should coincide with the oscillating term in the 
case of adiabatic transitions. The expression for the latter in terms of the functions e(t) and 
e'(t) can be easily derived from (1). It has the form of the term P£ sc (t,t ), eq. (24d), which 
allows to determine $22 (Mo) when P' = 0: 

<°(t, to) = f ^e 2 (f) + e«(f)df. (25) 

J to v 

In the majority of cases explicit expressions for $12 (Mo) anc ^ ^22 (Mo) do not exist. 
Nevertheless, eqs. (22) - (24d) permit to determine, in particular, some of the conditions 
under which the oscillating terms give negligible contributions in the probabilities P(v a ^ — > 
v a {p)) and P(v a (p) - > v p{a))\ they also allow one to obtain upper limits on these contributions 
when the latter are expected to be non- negligible. This can be done if the probability P' 
and the values of 9 m (t ) and 9 m (t), i.e., the values of the e(t) and e'(t) in the initial and 
final points of neutrino trajectory, are known. In this case the amplitude of the oscillations 
described by each of the terms given in eqs. (24a) - (24d) can be calculated. 

In many cases of physical interest one can use the Landau-Zener expression for the "jump" 
probability [|T2| , [T6|1 (see also |15| ): 



P'lz = e- 2 ™ , (26) 

where 

4n = An{t = t res ) = 2 f,fj. tr "| , (27) 

|^V 'res ) I 



9 There are several misprints in ref. (lj| (a list was given in ref. |p7| ). The three most relevant are: 
i) the overall minus sign in the right-hand side of eq. (40) should be canceled, ii) cos29 m (to) in 
eq. (41) should read sm26 m (to), and iii) a tx should be added in the right-hand side of the relation 
between the $i2(Mo) and the phase of the V\(to) — ► v^iti) transition amplitude, given in footnote 
9 ( c &i2(i)to) defined in |1^] is expressed in terms of arg (— A^) rather than of arg (A^)). Note 



that the probability amplitudes defined in [19] differ from those considered here by the phase factor 
exp (iLe(t')dt'), which is also reflected in the relation between $i2(Mo) and arg (A^) derived 
in PI. 



10 To our knowledge, the case of MSW transitions of solar neutrinos in the Sun (exponentially vary- 
ing density) is the only nontrivial and physically relevant one for which explicit analytic expressions 
for the phase functions $i2(Mo) an d < &22(Mo) were obtained so far. 



is the adiabaticity parameter, which is just the value of the adiabaticity function 

, s E™-E" 1 (e 2 (t) + e ,2 (t))l , N 

4n(t) = -?-. l - = 2 , ,\ \\ ;'\ - 28 

at the resonance point. In order for a given type of transition to be adiabatic, the inequality 
An(t) ^> 1 should be fulfilled at each point of the neutrino trajectory. In certain specific 
cases, such as MSW transitions of solar neutrinos in the Sun, the adiabaticity condition 
4n(t) ^> 1 is always satisfied if it is valid at the resonance point [0, i.e., if the inequality 
4n > 1 holds. 

The Landau-Zener expression for P', eq. (26), was derived assuming that e(t) decreases 
linearly with time (distance) while e'(t) is constant, on the neutrino trajectory. The limits 
of its applicability for describing the matter-enhanced neutrino transitions in a medium are 
well-known P, p!3|JT8| . It reproduces the "jump" probability rather accurately if the resonance 
region of the transition [[]] is sufficiently narrow, so that e'(t) does not change substantially 
and the change of e(t) in it can be well approximated by a linear function. Expression 
(26) cannot be used for description of nonadiabatic transitions when the point of neutrino 
production is located in the resonance region. 

In certain specific cases the analytic expression for the "jump" probability, derived for 
matter (electron, neutron number) density changing exponentially along the neutrino tra- 
jectory |I3|], P' exp i provides a more accurate description of the matter- enhanced neutrino 
transitions in a medium than P' LZ . For MSW neutrino transitions, eq. (6), the exponential 
density result for P' reads : 

„-27rn o (l-tan 2 0) „-27rno(tan~ 2 0-tan 2 6) a -2wr ^^- sin 2 9 „-27rr 4n^ 

™P = l_ e -2™ o (tan- 2 0-tan 2 0) = J _ ^2^4^ ' ^ ' 

Here r = \N(t = t Tes )/N{t = t res )\, N(t) = f t N(t) and N(t) = N e (t), i.e., r is the scale- 
height of the change of N e (t) in the case of interest. The analogous expression for the "jump" 
probability corresponding to RSFP in a constant magnetic field, P' exp {RSFP), can formally 
be obtained from eq. (29) by replacing Am 2 and 9 by Am 2 / cos 25 and 5 respectively, where 
tan 26 = 2/j, u B±/(Am 2 /2E), and by choosing the appropriate N(t) in the expression for r$ 
(in the case of eq. (7), for example, N(t) = N e (t) — N n (t)). 

Unlike expression (26), the one given by eq. (29) describes correctly, for instance, the 
(strongly) nonadiabatic MSW transitions of solar neutrinos in the Sun for values of sin 2 26 > 
(0.2 — 0.3) |J,13,18|. In the case of RSFP in a magnetic field which varies along the neutrino 



trajectory, one should use the value of the field at the resonance point, B±(t = t res ), in the 
expression for P^ xp (RSFP), i.e., in the definition of tan 25 given above. To our knowledge, 
the accuracy of the descriptions of the RSFP of solar and supernova neutrinos, based on the 
"jump" probabilities P' LZ and P^ xp (RSFP), has never been thoroughly tested for magnetic 
fields which change along the neutrino trajectory. 



ii 



This is the region of the neutrino trajectory around the resonance point, where sin 2 29 m (t) > 1/2. 



Let us note that independently of the value of P', the contribution of the oscillating terms 
P 1 °| c 34 (t, t ) in the probabilities P(v a (p) —* u a(/3)) an d P(v a (0) — > ^/3(a)) will be suppressed, 
as it follows from eqs. (24a) - (24d), if, for instance, sm29(t ) = (e.g., 9 (to) = tt/2) and 
sin29(t) = (e.g., 9(t) = 0). In the case of adiabatic transitions we have P' = and 
P 1 °^ 3 (t,to) — 0, while the oscillating term P£ sc (t,t ) is suppressed provided sin2#(£ ) — or 
sin29(t) =0. 

It is quite well-known (see, e.g., ||18|| ) that in the limits of sufficiently large and sufficiently 
small Am 2 / (2E) in the MSW case, the solar neutrinos oscillate in the Sun as in vacuum. It 
is instructive to see how expression (22) for P(y e — > v e ) reduces to the vacuum oscillation 
one in these two limits. For the solar neutrino MSW transitions in the Sun one has 9 m (i) = 9 
(see eqs. (6) and (9)) since N e (t) = at the surface of the Sun. If Am 2 /(2E) is sufficiently 
large, we have 4n > 1 and N r e es = Am 2 cos 29/ (2E^2G F ) > N e (t ) > N e (t'), t < t' < 
t, where N e (to) is the electron number density in the point of v e production in the Sun. 
Consequently, as it follows from eqs. (6), (9), (25) and (26) (or (29)), 9 m (t ) = 0, P' S 0, 
and $22 — Am 2 (t — t )/(2E). Using the above equalities and eqs. (22) - (24d) it is easy to 
convince oneself that P(v e —>■ v e ) indeed reduces to the vacuum oscillation probability. We 
get the same result when Am 2 /(2E) is sufficiently small, so that cos 29 m (t ) = —1 (-P 2 °| c 4 — 0) 
and 4no <C 1, i.e., the solar neutrinos undergo extremely nonadiabatic transitions. In this 
case eq. (29) implies []?' = cos 2 # and we have QHJ $i 2 + $ 22 = vr + Am 2 (t - t )/(2E). 
Therefore P(v e — > z/ e ) and P° sc reduce to the average probability and the oscillating term in 
the solar neutrino vacuum oscillation probability (see, e.g., p4f ): P(v e -^ z/ e ) = 1 — 1/2 sin 2 29 



and P° sc = l/2sm 2 29cos(Am 2 (t - t )/(2E)). Note that the nonadiabatic oscillating term 
P° sc converges to the vacuum oscillating term in this case. 

We would like to conclude this Section with the following remark. In certain cases of 
neutrino transitions in a medium one is faced with the problem of calculating the real part of 
the product of the two amplitudes A*(v a — > v a ) and A(v a —> z/g), arising as an interference 
term in the corresponding transition probability Q. Using eqs. (13a), (13b) and (19) - (21) 
it is not difficult to find the expression for the indicated product of amplitudes in terms of 
P', $i2 (Mo) and$ 22 (Mo): 

Re [A*(v a -> v a )A(v a -> v p )\ = ~ i\~ H ™s29 m (t ) sin29 m (t) + y/P'(l - P') cos29 m (t ) 



x cos 29 m (t) cos(<l>i 2 + $ 22 ) + y/P'(l - P') sin 29 m (t ) sin 29 m (t) cos($i 2 - $ 22 ) 

- -P'sin26' m (t )cos26' m (t) (cos2$ 12 + cos2$ 22 ) + -sm29 m (t )cos29 m (t) cos2$ 22 . (30) 

Obviously, the terms which depend on the phase functions $i 2 (t,t ) and/or $ 22 (t, t ) are 
oscillating terms. Note that the factors multiplying a given oscillating function (cos( < I ) i 2 + 



12 The Landau-Zener expression for P', eq. (26), leads to an incorrect result since it is not valid, 
in particular, for relatively small Am 2 /(2E) P Jl^Jlll . 

13 This term appears, for example, in the probability of combined MSW transitions and long wave 



length vacuum oscillations of solar neutrinos in the case of three flavour-neutrino mixing [ 28 1 . 
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$22), etc.) in the expressions for P(u a — > z/ a ;£,£ ) and Re [A*(u a — > u a )A(u a — > up)], i.e., 
the oscillation amplitudes, are different. The result (30) was used in the recent study [^] of 
the hybrid MSW + vacuum oscillation solution of the solar neutrino problem. It should be 
clear that using eqs. (13a) - (13d), (19) and (20) we can express any product of two of the 
amplitudes A^*'(u a — > v a (/3)) and A^*\vp —* vp(a)) in terms of P', $12, $22, m (to) and 9 m (t). 

3. MSW Transitions of Solar Neutrinos in the Sun and the Hydrogen Atom 

In the present Section we demonstrate that the second order differential equation for the 
amplitude A[y e — > v e ), which describes the MSW transitions of solar neutrinos in the Sun, 
coincides in form in the case of N e (t) changing exponentially along the neutrino path, with 
the Schroedinger equation for the radial part of the non-relativistic wave function of the 
hydrogen atom and we comment briefly on this interesting coincidence. 

Assuming that i) u a = u e , vp = z^( T ), ii) the initial conditions (A) are valid (A e (t,to) = 
A{y e — > Ue)), and substituting e(t) and e'{t) in eqs. (16) and (17) with their expressions 
given in eq. (6), we obtain from (17) a second order differential equation for the amplitude 
A[y e —>■ u e ), describing the MSW u e — > v^tr) transitions. According to the contemporary 



solar models |17 ], N e (t) changes approximately exponentially along the trajectory of a solar 



neutrino moving radially towards the surface of the Sun: 

iV e (£)=JV e (£ )expj-^j, (31) 

where (t — to) = d is the distance traveled by the neutrino in the Sun, N e (to) and r have 
been defined earlier, r ~ 0.1i? Q , R@ = 6.96 x 10 5 km being the solar radius. Introducing 
the dimensionless variable 

Z = ir V2G F N e (t )e~ t -^T, Z = Z(t = t ), (32) 

and making the substitution 

A e (t,t ) = A(u e - u e ) = (Z/Z ) c - a e-V-^fcW A' e (t,t ), (33) 



we find that the amplitude A' e (t, to) satisfies |]T3|,|T3|,|T9| the confluent hypergeometric equation 

{ Z^- 2 + (c - Z) ± - a\ A' e (t, to) = 0, (34) 



where []TQ | 



Am 2 2 Am 2 

a = 1 + 277) sin 9, c = 1 + ir Q . (35) 

u 2E 2E y J 

The equation (34) coincides in form with the Schroedinger (energy eigenvalue) equation 

obeyed by the radial part, ipki(f), of the non-relativistic wave function of the hydrogen 

atom |30|], \l/(r) = -ipki{ r )Yim{9' , 4>'), where r, 9' and <$' are the spherical coordinates of the 

electron in the proton's rest frame, I and m are the orbital momentum quantum numbers 

(m = —I, ■■■,1), k is the quantum number labeling (together with I) the electron energy 0, 



14 The principal quantum number is equal to (k + /) |?0| 
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Eki (Em < 0), and Yi m (9', </>') are the spherical harmonics. To be more precise, the function 
ip' kl (Z) = Z~ c l 2 e z l 2 ipki( r ) satisfies equation (34), where the variable Z and the parameters 
a and c are in this case related to the physical quantities characterizing the hydrogen atom: 



r 



Z = 2—^-E kl /Ej, a = a kl = l + l-^-Ej/E kl , c = q = 2(/ + l), (36) 

Go 

where ao = h/(m e e 2 ) is the Bohr radius and Ej = m e e 4 /(2/i 2 ) = 13.6 eV is the ionization 
energy of the hydrogen atom. It is remarkable that the behavior of such different physical 
systems as solar neutrinos undergoing MSW transitions in the Sun and the non-relativistic 
hydrogen atom are governed by one and the same differential equation. 

The properties of the linearly independent solutions of eq. (34), i.e., of the confluent 
hypergeometric functions, §(a,c; Z), as well as their asymptotic series expansions, are well- 
known p9f . Any solution of (34) can be expressed as a linear combination of two linearly 



independent solutions of (34), $(a, c; Z) and Z 1_c $(a— c+1, 2— c; Z), which are distinguished 
from other sets of linearly independent confluent hypergeometric functions by their behavior 
when Z — > 0: $(a', d\ Z = 0) = 1, a', d ^ 0,-1, —2, ..., a' and d being arbitrary parameters. 
Explicit expressions for the probability amplitudes A{y e — » u e ) and A{y e — + z^( T )) in terms 
of the functions $(a, c; Z) and $(a — c + 1, 2 — c; Z) were derived in [|19"1J2~61 . In the case of 
MSW transitions of solar neutrinos (N e (t) = 0) these expressions have an especially simple 
form: they are given by the corresponding vacuum oscillation amplitudes "distorted" by the 
values of the functions $(a', d\ Z) in the initial point of the neutrino trajectory, 

A(u e ^u li{T) ) = ± sin 20 j$(a-c,2- C ;Z )- e ^- t0 ^ $(a - 1, c; Z ) j, (37) 

etc., where Z , a and c are defined in eqs. (32) and (35). In the limit \Z \ — >• 0, which cor- 
responds to zero electron number density, expression (37) reduces to the one for oscillations 
in vacuum. 

It is well-known that the requirement of a correct asymptotic behavior of the wave func- 
tion ijjki{ r ) at large r leads to the quantization condition for the energy of the electron, E k i, 
in the hydrogen atom 0: E M = -E 1 j{k + I) 2 , (k + I) = 1, 2, ... (I = 0, 1, 2, ..., (k + I) - 1). 
Technically, the condition is derived by using the asymptotic series expansion of the con- 



fluent hypergeometric functions in powers of the argument Z [2E] (one has Z — > oo when 
r — > oo, see eq. (36)). The same asymptotic series expansion in the case of the solutions 
describing the MSW transitions of solar neutrinos in the Sun (we have \Zq\ ^ 520 in this 
case ||19| ) permitted to obtain expression (29) for the "jump" probability P' [I3|, as well as 
explicit expressions for the phase functions $12(^,^0) and $22 (Mo) EH- 

4. Conclusions 

We have derived an exact universal analytic expression for the probability of two-neutrino 
matter-enhanced transitions in a medium (MSW, RSFP, induced by neutrino FCNC inter- 
action, etc.). The probability is expressed in terms of three real functions of the parameters 
characterizing the neutrino transitions: the "jump" probability P' and two phases (angles) 
$12 and $22- The latter are responsible for the oscillatory dependence of the probability on 
the parameters of the problem being investigated (neutrino energy, matter density, magnetic 
field strength, etc.). Although in the majority of cases of physical interest the two phase 
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functions $12 and $22 are not known, the amplitudes of the oscillations due to $12 and $22 
are functions only of the "jump" probability and the values of the neutrino mixing angle in 
matter at the initial and final points of the neutrino trajectory, i.e., of quantities which are 
typically estimated or calculated in the studies of the neutrino matter-enhanced transitions. 
Our results can be used, in particular, in the investigations of matter-enhanced transitions 
of solar and supernova neutrinos. 

We have noted also that the second order differential equation for the MSW probability 
amplitude of solar v e survival in the Sun coincides in form, for solar matter (electron number) 
density changing exponentially along the neutrino trajectory, with the Schroedinger equation 
for the radial part of the non-relativistic wave function of the hydrogen atom. The equation 
is of a confluent hypergeometric type and the asymptotic series expansion of its solutions 
plays important role in the description of both physical systems: it permits to obtain, for 
instance, the quantization condition for the electron energy levels in the hydrogen atom and 
the analog of the Landau-Zener "jump" probability for exponentially varying density in the 
case of MSW transitions of solar neutrinos in the Sun. 
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